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1 Introduction 

Several theoretical difficulties of the standard big bang cosmology such as the horizon and flatness 
problems can be elegantly solved by inflation [1]. In fact, the existence of the inflationary era in 
the early Universe is strongly supported by the observations [2] ; the density perturbations extending 
beyond the horizon at the last scattering surface can be interpreted as the evidence for the accelerated 
expansion in the past. 

The study of density perturbations such as isocurvature perturbations, non-Gaussianity, tensor- 
mode, and their effects on the cosmic microwave background (CMB) power spectrum is a powerful 
diagnostic of the mechanism that laid down the primordial density fluctuations, but it is not enough 
at present to pin down the model. This is partly because of our ignorance of thermal history of 
the Universe beyond the standard big bang cosmology, especially concerning how the Universe was 
reheated. 

Whereas one of the plausible explanation for the density perturbation is the quantum fluctuations 
of the inflaton from the minimalistic point of view, it may be that there are many other light scalars in 
nature, one of which is responsible for the observed density perturbation via the curvaton [3, 4, 5, 6] 
(or its variant, e.g. modulated reheating [7, 8]) mechanism. In fact, there are many moduli fields 
that necessarily appear at low energies through compactifications in the string theory. Most of them 
must be stabilized in order to have a sensible low-energy theory, but some of them may remain 
relatively light, and therefore are a candidate for the curvaton. Interestingly, there is an argument 
that the string theory contains a plenitude of axions, the so called "string axiverse." [9] We shall see 
later that the axion is indeed a plausible candidate for the successful curvaton.^ 

* Throughout this paper the axions refer to imaginary components of moduH fields, which have a sinusoidal potential. 
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One of the distinguishing features of the curvaton mechanism is that it can generate the density 
perturbation with large non-Gaussianity. If any primordial non-Gaussianity is found by the Planck 
satellite, it would exclude a simple class of inflation models as the origin of the entire density 
perturbation, and therefore, it has a tremendous impact on our understanding of the early Universe. 

Recently, the present authors studied non-Gaussianity generated by the curvaton mechanism in a 
great detail, and developed a formalism to calculate the density perturbation for a generic curvaton 
potential [10]. We pointed out that the curvaton should be located at a potential with negative 
curvature during inflation, and in particular it must be close to the local maximum ("hilltop") 
of the potential, in order to generate a red-tilted density perturbation spectrum which is strongly 
favored by the recent observations [2]. Interestingly, we found that, even if the curvaton dominates 
the Universe, the non-Gaussianity parameter /nl is positive and gets enhanced logarithmically in 
the hilltop limit, and therefore /nl of 0(10) is a robust prediction of the hilltop curvaton. Applying 
our formalism to the axionic (or pseudo-Nambu-Goldstone) curvaton with a sinusoidal potential, 
we found that /nl can be as large as about 30, which is realized for the curvaton mass of order 10 
TeV and the decay constant of order the GUT scale. In this analysis we fixed the scalar spectral 
index Ug = 0.96 for simplicity. The mild increase of the non-Gaussianity in the hilltop limit is 
originated from the fact that the density perturbation generated by the curvaton is enhanced. This 
enhancement is due to non-uniform onset of curvaton oscillations [11, 10]. This result should be 
contrasted to a simple curvaton model with a quadratic potential, which predicts a negative /nl of 
order unity in the case that it dominates the Universe. 

In this paper, we extend our previous work on the non-Gaussianity generated by the axionic 
curvaton with the hilltop initial condition. We will discuss its dependence on the scalar spectral 
index, and also scan the curvaton parameters, namely, the mass and the decay constant. Interest- 
ingly, we find that /nl is bound as /nl ^ 30 for = 0.94-0.99, and the maximal non-Gaussianity 
is realized for the curvaton mass 10 TeV and the decay constant of order the GUT scale. (If reheat- 
ing happens prior to the curvaton oscillation, then the bound becomes /nl ^ 40.) Furthermore, 
/nl = 20-40 is realized for a wide range of parameters, the curvaton mass rn-o- = 10- 10^ GeV and 
the decay constant / = 10^^ - 10^ GeV. One of the plausible candidates for such an axionic curvaton 
is an imaginary component of the moduli (i.e., axions) with mass of order 10 -100 TeV and decay 
constant of 10^^"^ ''GeV. The moduli fields are stabilized by the non-perturbative effect and the 
supersymmetry (SUSY) breaking, and it is plausible that the moduli mass is closely related to the 
SUSY breaking scale in the visible sector. Intriguingly, such several tens TeV SUSY breaking scale 
is consistent with the recently discovered Higgs boson mass of 125- 126 GeV [12, 13]. 

The rest of the paper is organized as follows. After briefly reviewing density perturbations from 
general curvatons in Section 2, then in Section 3 we discuss axionic curvatons in detail. We then 
give discussions and conclusions in Section 4 and 5, respectively. 

The appendix discusses an extreme case where the curvaton drives a second inflationary period. 
After analytically computing density perturbations from inflating curvatons in general, we then 
apply the discussions to axionic curvatons. We find that the non-Gaussianity turns out to be rather 
small when the axionic curvaton drives a second inflation. 



2 



2 Review of Curvatons with a General Potential 



In the curvaton mechanism, the hght curvaton field acquires super-horizon field fluctuations during 
inflation. The density perturbations are produced in the post-inflationary era, as the curvaton 
oscillates and its energy density relatively grows compared to other radiation components. In this 
section we give a brief review of density perturbations generated by a curvaton a with a generic 
effective potential V{(t). We refer the reader to [10] for detailed derivation of the following results. 

2.1 Density Perturbations from Curvatons 

The density perturbations generated by curvatons depend on the curvaton dynamics during and 
after inflation. In the simple curvaton model with a quadratic potential, the curvaton dynamics is 
determined by the curvaton mass and the initial deviation from the origin. If the mass is much 
smaller than the Hubble parameter during inflation, the curvaton does not evolve until it starts to 
oscillate, and the resultant density perturbation is given in a rather simple form. However this is no 
longer the case for a general curvaton potential. In particular, as we shall see shortly, the curvature 
of the potential should be negative and non-negligible in order to account for a red-tilted density 
perturbation spectrum, and therefore, the curvaton significantly evolves after inflation, affecting the 
density perturbation. 

If the curvaton potential V{a) has no explicit dependence on time, then the curvaton dynamics 
prior to the oscillation can be tracked by the attractor solution 



cHa = - V, with 



3 (during inflation with H ~ const.) 
9/2 (matter domination) (2-1) 
5 (radiation domination) 



which is a good approximation while \ V" /cH'^\ <^ 1. Here, a prime denotes a derivative with respect 
to a, an overdot a time derivative, and H = a/ a. Setting the minimum of the potential about which 
the curvaton oscillates to cr = 0, the onset of the oscillation can be defined as when the time scale 
of the curvaton rolling becomes comparable to the Hubble time, i.e. 



(T 



1. (2.2) 



Ha 

Then the Hubble parameter at the time is obtained as 

hL = (2.3) 

where the subscript "osc" denotes values at the onset of the curvaton oscillation, and c is a constant 
depending on whether reheating (= inflaton decay, at treh) is earlier /later than the onset of the 
curvaton oscillation (corresponding to c in the attractor (2.1) right before the oscillation): 



9/2 (treh>to 
5 {treh < 



(2.4) 



Let us here summarize simplifying assumptions concerning the evolution of the energy densities 
of the curvaton and the inflaton. We assume the curvaton potential to be well approximated by 
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a quadratic one around its minimum so that the curvaton oscillations are sinusoidal.^ Then the 
curvaton energy density redshifts similarly to nonrelativistic matter after the onset of the oscillations 
until the curvaton decays into radiation. On the other hand, we consider the inflaton to behave as 
matter from the end of inflation until reheating when it decays into radiation. The energy density 
of the curvaton before the beginning of its oscillation is assumed to be negligibly tiny compared to 
the total energy of the Universe, having little effect on the expansion history. 

Supposing the curvaton field fluctuations to be nearly Gaussian with Vsa{k) = {H\k=aH f'^'^)'^ at 
the time when the comoving wave mode k exits the horizon, then using the (5AA-formalism [14, 15, 
16, 17], the power spectrum of the density perturbations at the CMB scale is expressed as [10] 



with 

(i-xw^yr^i^-'-^^]^. (2.6) 



(2.7) 

dec 



da^ 4 + 3r i ^(o-qsc) o-qsc j V'{a^) 

Here, the subscript * denotes values when the CMB scale exits the horizon, and r is the energy 
density ratio between the curvaton and radiation (which originates from the inflaton) upon curvaton 
decay 

r = — 

Pr 

The function X denotes effects due to the non-uniform onset of the curvaton oscillations (which are 
absent for a purely quadratic curvaton potential), defined as follows: 

^Wosc) = —f ^ — r^TT ^ 1 , (2.8) 

2(c - 3) V V'{aosc) ) 

where the constant c is given in (2.4). 

Prom the above expressions, the spectral index of the linear order perturbations follows as 

The recent observations strongly suggest that the density perturbation power spectrum is red-tilted, 
ris = 0.968±0.012 [2]. As we shall see later, this requires that the curvaton potential be tachyonic and 
the size of the curvature must be of order 10 % of the Hubble parameter during inflation, unless the 
inflaton is allowed to take super-Planckian field values, or some special configurations are arranged 
in the inflationary setup. 

Curvatons also generate local-type^ bispectrum, whose amplitude is represented by the non- 



^Cases with non-sinusoidal oscillations are discussed in Appendix B of [10]. 

^Strictly speaking, bispectra from curvatons have shapes similar to, but may not exactly be of the "local form" [18], 
especially when /nl is strongly scale-dependent [19, 20, 21]. However we note that for axionic curvatons with sinusoidal 
potentials, the running of /nl is tied to the running of the spectral index, and thus strictly constrained to be small by 
current observations [211. 



4 



linearity parameter /nl. This is given by 

_ 40(1 + r) ^ 5(4 + 3r) [ V'{aosc) 3X(aosc) 



3r(4 + 3r) 



6r 



V'(crosc) 



(l-X(cTosc))"'X'(fT, 



y^(aosc) 3X(aosc) r' / y^^(cTosc) y^(aosc)' 3X^(aosc) ^ 3X(aosc) 



osc 7 



(To 



(2.10) 



+ T777 T - (1 - ^(i7osc))t777 



l^'((Tosc) 



l^'((Tosc) 



In the limit of X((Tosc) — ^ 0, this reproduces the known result for the curvaton with a quadratic 
potential. 

Let us also rewrite the energy density ratio r (2.7) in terms of the inflaton and curvaton param- 
eters: 

„ 4/3] 



Max. 



X Min. 1, 



r 



1/2 



Hole j [ 'iM^Holc^a' 



/2 



X Min. 1, 



0-1/2 



2.11 



where Mp ~ 



2.4 X 10^® GeV is the reduced Planck mass, and the first and second terms in the Max. 
parentheses correspond to the curvaton being subdominant and dominant at its decay, respectively, 
while the Min. parentheses are due to whether the onset of oscillation is after or before reheating. 
and To- are constants that denote, respectively, the decay rates of the inflaton and the curvaton. We 
note that in obtaining the above results, we have adopted the sudden decay approximation where 
the scalar fields suddenly decay into radiation when H = T. 

Finally, the curvaton field value at the onset of the oscillations cJosc is obtained by integrat- 
ing (2.1) 

2.12 



O-osc 

V' 



ml, 2c{c-mL 



solving which gives cTosc as a function of a^J Here, A/"* is the number of e-folds during inflation 
between the horizon exit of the CMB scale and the end of inflation, c is given in (2.4), and Hi^f is 
the inflationary Hubble scale (we are assuming a nearly constant Hubble parameter during inflation, 
thus -ffinf ^ H^:). 

Therefore by combining the above expressions, one can compute the density perturbations from 
a curvaton with a generic potential V{(j), given the curvaton field value at the CMB scale horizon 
exit the decay rates of the inflaton and curvaton T^, the inflationary scale -ffmf, and the 
duration of inflation A/"* . 



2.2 Case Study: Hilltop Curvatons 

As an example that will be relevant for analyzing axionic curvatons in the next section, here let us 
apply the above generic results to a curvaton located at the hilltop, whose potential around (Tosc and 

^When (2.12) admits as solutions for (Josc both positive and negative values, one should take the sign of aosc to 
match with that of a,. 



5 



£7=1, is well approximated by 

V{o) = Vo-\m\a-a^)\ (2.13) 

where m, (Tq, and Vo(> 0) are constants. Without loss of generality, we assume < o"osc < cr* < ctq- 
Then one can check that when the curvaton is close enough to the hilltop to satisfy 

r, (2.14) 

then the resulting power spectrum (2.5) and the non-Gaussianity (2.10) take the form 

1/2^ 3r_ao^ao^^^ (2.15) 



^ 4 + 3r (To — cT* 27r(T, 



OSC 



with spectral index (2.9) 



^ 5(4 + 3r) dose 

JNL ^ 7^ , (2.1bj 

18r do - o-Qsc 



i?* 2m2 ^ ^ 



The equation (2.12) which relates o"* and o"osc gives 

Inf , ^ , , (2.18) 

\(To-CJosc/ 2(c - 3) do - O"osc 

where we dropped the H?^^ contribution in the right hand side from (2.14) and by assuming 
™^ l^ini i$ 10^^- initial value is shifted towards the hilltop, fiosc approaches do much 

slower than cr=K does since the left hand side is logarithmic. Therefore as one approaches the hilltop, 
Vc^ (2.15) blows up due to the enhancement factor (do — dosc)/(co — c*)5 while /nl (2.16) increases 
slowly. We also note that the value of /nl is greater than one even when r ^ \, from (2.14). The 
extreme amplification of the linear perturbations corresponds to the curvaton taking longer time to 
start its oscillation when starting closer to the hilltop. 

Before ending this section, we should remark that in the extreme hilltop limit, the approxima- 
tion (2.1) for the curvaton dynamics mildly breaks down before the curvaton starts to oscillate. This 
gives rise to errors of 0{1) for the above results in this limit. However, the above analytic expres- 
sions suffice for our order of magnitude estimations on axionic curvatons in the next section. We 
will also carry out numerical computations when further accuracy is required, e.g., when calculating 
predictions on /nl- 



3 Axionic Curvatons 

Now let us move on to the investigation of axionic curvatons, which is the main topic of this paper. 
As was explained in the introduction, we focus on the case where the curvaton is a pseudo-Nambu- 
Goldstone boson of a broken U(l) symmetry, possessing a periodic potential of the form 



V{a) = 



1 



cos 



d 



(3.1) 
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where / and A are mass scales. Without loss of generality, we restrict the initial field value to lie 
within the range < o"* < vr/. The curvaton's effective mass at the potential minimum is denoted 
by 

= —. (3.2) 

Then supposing that the coupling of the axionic curvaton with its decay product is suppressed by 
the symmetry breaking scale /, the curvaton decay rate takes the value 

where the constant /3 is naively of order unity. In the following, we ignore the time-variation of 
the Hubble parameter during inflation, and especially, neglect the H contribution to the spectral 
index (2.9). In other words, we do not consider inflationary models with rather large \H /H'^\ which 
requires super-Planckian field ranges'^ or some special configurations. Hence the axionic curvaton 
need to be located beyond the inflection point during inflation, i.e. 0.5 < cr*/vr/ < 1, in order to 
source a red-tilted power spectrum. 

The axionic curvaton with cr* <^ nf whose potential is well approximated by a quadratic was 
studied in [23], and the whole potential including the hilltop region was investigated in [10]. There 
it was shown along the line of discussion in Section 2.2, that unless the axionic curvaton is initially 
located close to the hilltop, both the inflation and reheating scales need to be very high. For e.g., 
for (T^/tt/ = 0.75 to satisfy both the WMAP normalization ~ 2.42 x 10~^ and the spectral 
index Ug ~ 0.96, then Hi^{ > 10^^ GeV and p^^^ > 10^^ GeV are required, where preh represents the 
radiation energy density at the reheating. This is because the spectral index of order 1 — ~ 0.01 
requires a rather large curvaton mass m„ ~ 0.1-ffinf, forcing the curvaton to start its oscillation soon 
after the end of inflation. Hence without high inflation and reheating scales, the curvaton cannot 
even come close to dominating the Universe to source measurable density perturbations.^ The story 
is quite different for an axionic curvaton in the hilltop region, where the onset of the oscillation is 
delayed and curvaton domination is allowed with lower inflation/reheating scales. This, together 
with the amplification of the linear perturbations in the hilltop limit (cf. discussions around (2.18)), 
makes axionic curvatons compatible with many orders of magnitude of the inflation and reheating 
scales. 

In light of the above considerations, in this section we elaborate on axionic curvatons in the hilltop 
region, which dominate the Universe before decaying into radiation. We will find that this particular 
limit of axionic curvatons has interesting predictions, especially in terms of the non-Gaussianity. 



Assuming single-field canonical slow-roll inflation, the Lyth bound [22] relates the time-variation of the Hubble 
parameter with the inflaton field (j) range as 

where Mp is the reduced Planck mass and JV the e- folding number. Thus \ H/H^\ as large as to give sizable contribution 
to the spectral index (2.9) whose typical value is ~ 0.968 (WMAP central value) normally requires a super-Planckian 
field range for the infiaton. 

^The curvaton's effective mass during infiation is decoupled from the mass at the potential minimum (3.2) when 
the curvaton is close to the inflection point, i.e. at/vr/ ~ 0.5, however in such case even higher inflation/reheating 
scales are required. 
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3.1 Parameter Space in the Hilltop Regime 

The axionic curvaton model has five free parameters, which are the symmetry breaking scale /, the 
effective mass mo- = A^//, the curvaton field value at CMB scale horizon exit cr*, the inflationary 
scale -ffinf, and the inflaton decay rate T^- However, since we are focusing on a curvaton that 
dominates the Universe before it decays, as long as there exists a parameter window which allows 
r ^ 1, the cosmological observables do not depend on the explicit value of r or Fq-. In this sense, 
the dominant axionic curvaton is actually a four parameter model. 

Strictly speaking, there are three more parameters: the e-folding number A/"* between the CMB 
scale horizon exit and the end of inflation, the constant c (2.4) representing whether tj-^h ^ iosc 
(though this is determined when the other parameters such as F,^ are fully given), and /3 in (3.3) 
parameterizing the curvaton decay rate. A/"* determines how much the curvaton rolls during inflation 
(cf. (2.12)), however such rolling is negligible compared to that in the post-inflationary era as seen 
in (2.18), and thus has little effects on the model. Hence we simply flx the e-folding number to 
A/"* = 50 in the following discussions. As for c, whether reheating happens before/after the onset of 
the curvaton oscillations do not affect the allowed parameter window for / and ttt-o-, but give slightly 
different predictions on /nl- This will be discussed in Section 3.2. The parameter /3 for the decay 
rate is set to unity in the following, and implications of (3 taking other values are also discussed 
later. 

Out of the four parameters, ifjnf and cr*// can be fixed from the WMAP normalization 

PS 2.4 X 10"^ (3.4) 

and also from requiring the spectral index to be consistent with the WMAP bound 

Us ^ 0.96. (3.5) 

Later on we will see that the detailed value of the spectral index, as long it is not so close to unity, 
only have minor effects on axionic curvatons. Hence we are left with two parameters for the axionic 
curvaton / and ma- Order of magnitude constraints on these parameters can be obtained using the 
analytic formulae in Section 2 (or in Section 2.2), which we present in Figure 1. The yellow region 
corresponds to the allowed window for a dominant axionic curvaton in the hilltop, i^inf and u*// 
are fixed to appropriate values by the observational constraints (3.4) and (3.5) at each point in the 
window, as indicated in the upper figures showing their contour lines. This also fixes the curvaton 
decay rate via (3.3), cf. lower left figure. (The relativistic degrees of freedom is fixed to 5* = 100 
upon drawing the T^ec contours). On the other hand, the infiaton decay rate F^ is not fixed at each 
point but is allowed to take values within a certain range, as we will soon explain. The lower right 
figure shows contour lines for the non-Gaussianity /nL) which is typically a few tens. Here we note 
that since the analytic formulae in the previous section can contain 0(1) errors (cf. discussions at 
the end of Section 2.2), the /nl values have been computed numerically. We have only shown the 
/nl contours inside the allowed window for hilltop axionic curvatons, but the values can be modified 
at regions very close to the boundaries. 

We find that the allowed window is constrained by the following four conditions: The upper edge 



8 



(green line) is set by the requirement that the curvaton initiahy hes in the hihtop regime, 

> 0.9. (3.6) 

Recah that a non- hihtop axionic curvaton can work only with very high inflation/reheating scales. 
The right edge (blue line) denotes the requirement that the curvaton be subdominant until it starts 
its oscillation, 

F(cTosc) < 0.1 X 3M^Hi,. (3.7) 

When going beyond this boundary, the curvaton starts to drive a secondary inflation. A rather strict 
relationship between iJm and ex*// is required for such inflating curvatons to work, as is discussed 
in Appendix A. The lower edge (orange line) requires the curvaton to decay at temperatures higher 
than 5 MeV in order not to ruin Big Bang Nucleosynthesis (BBN) [24, 25, 26, 27], i.e. 

SM^r^ > ^g,(5MeV)^ (3.8) 

with the relativistic degrees of freedom g^, = 10.75. Finally, the left edge (red line) follows from the 
dominant condition 

r > 10. (3.9) 

These conditions give the most stringent constraints on the hilltop axionic curvaton, and other 
requirements for a consistent curvaton scenario are satisfled in the window bordered by (3.6) - (3.9). 

Let us also lay out such satisfled conditions: Firstly, the curvaton energy density is negligibly 
small during inflation, 

V{a,) < SM^Hl^. (3.10) 

Moreover, quantum fluctuations during inflation should not make the curvaton jump over its poten- 
tial minimum in order to avoid the resulting density perturbations from being highly non-Gaussian, 
or over the maximum to avoid domain walls. 

In the hilltop region, the classical rolling becomes suppressed, which can compete with the quan- 
tum fluctuations during inflation. The curvaton's classical rolling dominates over the quantum 
fluctuations if 



3 H? 



inf 



<. 1, (3.12) 



2iT V'(a^, 

where the curvaton is considered to slow- roll due to (3.10) and the lightness condition that follows 
from the spectral index (3.5). Furthermore, the curvaton decay should happen after reheating and 
the onset of the oscillations, 

Tcr < r,^, -ffosc- (3.13) 

^" If the curvaton is subdominant at the decay, namely, r ^ 1, the non-Gaussianity tends to be too large, and one 
has to tune the decay rates of the curvaton and the inflaton so that r > 0.01 is realized in order to be consistent with 
the observations. We have nothing new to add to this possibility in our context. 

^^The randomized case of axionic curvatons at the potential minimum is discussed in [23]. 
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The mass nia is required to be larger than the curvaton decay temperature, in order to avoid possible 
backreaction effects to the curvaton's perturbative decay (see e.g. [28, 29, 30]). Assuming instant 
thermalization, this condition is written roughly as 



As for the inflaton sector, reheating (= inflaton decay) should happen subsequent to inflation, while 
an upper bound on the inflationary energy scale is given by constraints on primordial gravitational 
waves. The 7-year WMAP+BAO+i?o gives Vt/'Pi^ < 0.24 (95% CL), which translates into^^ 



Let us repeat that all the requirements (3.10) - (3.15) are satisfied in the allowed window of Figure 1. 

We should also remark on the constraints on the reheating scale before ending this subsection. 
As we have noted above, dominant axionic curvatons are insensitive to the explicit value of T^. The 
only constraints on are that the inflaton should decay after the end of inflation (3.15) but before 
the curvaton decay (3.13). For the case of tosc < ^reh, the dominant condition (3.9) sets an additional 
lower bound on F,^, cf. (2.11).^'^ The inflaton decay rate should take values within these bounds at 
each point of the allowed window in Figure 1 . We note that the contour lines of various quantities in 
the figures are obtained assuming that the values of F,^ at each point do not saturate the lower / upper 
bounds set by the above requirements. If, for example, F^ takes lowest possible values saturating 
the dominant condition (3.9) for the tosc < ^reh case, then r becomes as small as ~ 10, slightly 
modifying /nl from the shown values. 

3.2 Dependence on Various Parameters 

The spectral index 7x^ — 1 = — 2m^/3i7^£ fixes the inflationary scale i^inf proportional to the curvaton 
mass rrio-, as is shown in the upper right figure. The rather wide range allowed for is translated 
into the axionic curvaton being compatible with inflationary scales with many orders of magnitude. 

The upper left figure shows the curvaton's initial value a^:, which comes closer to the hilltop vr/ 
towards the lower right corner of the allowed region. However when so close to the hilltop such that 
a^:/irf ^ 1 — 10^^^, then the axionic curvaton either ruins BBN, or drives a secondary inflation. 
For dominant curvatons r ^ 1, the non-Gaussianity /nl is determined by how close the curvaton 
initially is to the hilltop (cf. (2.16) and (2.18)), thus the /nl contours run parallel to those of 
(vr/ — (T*)/7r/. Here, recall that when shifting cr* towards the hilltop, dose increases much slower 
than fj* does. This leads to a mild increase of /nl, whose largest possible value is ~ 30 for the case 

^■^When the inflation scale is high enough to saturate the bound (3.15), then depending on the inflationary mechanism, 
one can expect to have a contribution to the spectral index from a non- vanishing H/H^ (see also Footnote 8), as well as 
the central value of the spectral index bounds (3.5) being slightly shifted. We ignore such effects for axionic curvatons, 
since they can modify the results only at the vicinity of the upper right corner of the allowed window in Figure 1. 

^■^For the case of tosc < ireh, the left red edge actually denotes where the upper bound on set by F^ < Hobc and 
the lower bound from the dominant condition (3.9) take the same values. In other words, tosc < ^rch and (3.9) are 
incompatible beyond the red line. On the other hand, for the tosc < troh case, r is independent of F0, which allows 
one to set the dominant condition (3.9) independently of However we note that (3.9) produces the left red edge 
at the same place on the / — rricr plane for both trch ^ tosc cases. 




(3.14) 



F^ < /?inf < 1.3 X lO^^GeV. 



(3.15) 
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Figure 1: Parameter space for a dominant axionic curvaton in the hilltop. The allowed window is 
shown as the yellow region, which is bordered by the hilltop condition (3.6) (upper green bound- 
ary), the requirement that the curvaton does not inflate the Universe (3.7) (right blue), BBN con- 
straint (3.8) (lower orange), and the dominant condition (3.9) (left red). These conditions have been 
adopted in order to study the peculiar behavior of the hilltop curvaton. However the curvaton mech- 
anism can still work when relaxing some of them, see the discussions in Appendix A and footnote 10. 
The contour lines on each figure denote the following quantities. Upper left: The curvaton value at 
CMB scale horizon exit (vr/ — o"*)/7r/. Upper right: Inflationary scale H[^f in units of GeV. Lower 
left: Decay temperature T^ec of the curvaton in units of GeV. Lower right: Non-Gaussianity /nl for 
the case of tosc < ^rch- The values of /nl slightly increases when tosc > ^reh- 
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of tosc < irchi and ~ 40 for tosc > ^rch- We also note that the regime beyond the right blue edge 
corresponds to axionic curvatons driving a second inflationary stage. However, such inflating axionic 
curvatons produce rather small non-Gaussianity, as discussed in detail in Appendix A. 
Now let us discuss the model dependence on other parameters. 

Reheating Before/During Curvaton Oscillations 

Whether reheating happens before or during the curvaton oscillation only slightly modify the cur- 
vaton velocity prior to the oscillation. This does not affect our order of magnitude estimation on 
the allowed window in the / — nia plane, except for that the case of trch < ^osc restricts the inflaton 
decay rate to lie within a rather narrow range [10]. The ct^k/tt/, -ffmf, and T^ec contours are also 
nearly the same for the two cases, however we note that /nl can be slightly larger when tj-eh < ^osc- 
This is because a radiation dominated Universe allows a to roll less compared to when dominated 
by matter, and thus slightly makes dose closer to the hilltop. The lower right figure shows the /nl 
contours for ireh > ^osc, but the case of treh < *osc increases the /nl values by up to ~ 10. 

Spectral Index 

The allowed window and non-Gaussianity are insensitive to the explicit value of Ug (whether 
is, say, 0.94 or 0.98). However, if the spectral index is as close to unity as Ug > 0.99, then the 
curvaton potential is required to be so flat such that the quantum fluctuations can dominate over 
the curvaton's classical rolling during inflation in the hilltop regime. When increasing Ug beyond 

0. 99 towards unity, the condition (3.12) is violated flrst in the lower right corner of the allowed 
window in the / - nia plane, and eventually in the entire window at Ug > 0.999. Normally the model 
loses precise predictions if the quantum fluctuations dominate over the classical rolling. However 
we expect that the predictions for hilltop curvatons are not affected much by such quantum jumps 
during inflation, since it is the non-uniform onset of the curvaton oscillation that mainly generates 
the linear and second order perturbations, and also because (3.11) is satisfled even for Ug w 0.999, 

1. e., the quantum jumps (even when they dominate the curvaton dynamics) do not drastically change 
the curvaton position during inflation. We leave this question for future work, and let us close this 
paragraph by stating that as long as Ug < 0.99, the detailed values of the spectral index has little 
effect on axionic curvatons. 

Curvaton Decay Rate 

We have been setting (3 as unity in the curvaton decay rate (3.3). A further suppressed To- (i.e. 
smaller /3) delays the curvaton decay, thus makes the BBN constraint (3.8) more stringent, while 
making it easier for the curvaton to dominate the Universe and relaxes the dominant condition (3.9). 
For example, /3 = 10~^ tightens the lower edge (orange boundary) of the allowed window in Figure 1 
by A(logio "^cr) ~ 1, but pushes out the left edge (red) slightly (i.e. does not change the order of /). 
It should be noted that the tightening of the BBN constraint results in decreasing the largest possible 
value for /nl, as can be seen in the lower figure. For /3 = 10~^, the maximum /nl is about 27. 
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4 Discussion 



The upshot of our results is that the non-Gaussianity parameter /nl is generically of 0(10) in the 
hihtop curvaton, even when the curvaton dominates the Universe. In particular, as one can see from 
Fig. 1, /nl = 20-40 is realized for the curvaton mass nia = 10-lO^GeV and the decay constant 
/ = 10^2. 1017 GeV. 

What is the plausible candidate for the axionic curvaton? Interestingly, there are many moduli 
fields (r) in the string theory, and they are massless at the perturbative level because of the shift 
symmetry, 

T T + ia, (4.1) 

where q is a real transformation parameter. After the moduli fields are stabilized by non-perturbative 
effects and SUSY breaking, the imaginary components of the moduli fields, namely the (string) 
axions, acquire a sinusoidal potential like Eq. (3.1). The symmetry breaking scale / is naively 
expected to be of order the GUT or Planck scale. Thus, the string axion is one of the plausible 
candidates for the axionic curvatons.^'^ 

Recently, the standard-model like Higgs boson was discovered by the ATLAS and CMS exper- 
iments [12, 13]. The observed Higgs boson mass is about 125-126GeV, which suggests relatively 
high-scale SUSY [31, 32]. Depending on the ratio of the up- and down-type Higgs boson VEVs, the 
SUSY breaking scale ranges from lOTeV up to several tens PeV. While the axion mass crucially 
depends on the stabilization mechanism, it is related to the gravitino mass in a KKLT-type stabi- 
lization [33], and so, it is conceivable that the axion mass is not many orders of magnitude different 
from the suggested SUSY breaking scale in the visible sector. It is intriguing that the axion with 
mass of this order can generate a large non-Gaussianity within the reach of the Planck satellite. 

The initial position of the curvaton must be very close to the hilltop of the potential. If some 
symmetries are restored at the maximum of the potential, the curvaton sits initially very close to the 
hilltop without any fine-tuning. This is possible if one considers a moduli space spanned by multiple 
scalar fields [34]. To be concrete, we consider a supersymmetric theory with the superpotential, 

W = 5(^2 _ ^2 _ ^2)^ (4 2) 

where S, x and (f) are chiral superfields, and ;U is a mass scale. In the supersymmetric limit, there is 
a moduli space is characterized by 

x' + 0' = m', (4.3) 

where it should be noted that both x a-nd <j) are complex scalar fields. The scalar potential vanishes 
in the moduli space. There are two special symmetry-enhanced points, i.e., x = and = 0. 
The degrees of freedom orthogonal to the moduli space are heavy, and can be integrated out. For 
instance, x is heavy at (/> ~ 0, one of the symmetry enhanced points, and we can erase x by using 
(4.3). In order to see that the potential has extrema at those symmetry enhanced points, let us 
introduce a soft SUSY breaking mass, m^J^J^^ which lifts the moduli space. A similar soft SUSY 

The real component may play a role of the inflation, in which case the moduli explains both the inflation and the 
origin of density perturbations. 
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breaking mass can be introduced for (p, but it does not change the argument. Since it is (p that is 
Hght at w 0, the effective potential can be written as 

Ves = mV'-y^'l, (4.4) 

where we have minimized the angular component of (j), and defined (p = Thus, </> = is the local 
maximum if > 0. Note that one should write the effective potential in terms of x at <^ ~ 
since (j) becomes heavy and it is x that is light. Then the potential is simply given by m^|xP, which 
clearly shows that the potential is minimized at x = (or if = 

Now let us discuss other cosmological issues. In order to have successful cosmology, it is neces- 
sary to generate a right amount of baryon asymmetry and dark matter. Since the baryonic/CDM 
isocurvature density perturbation is tightly constrained by observations, it also limits possible baryo- 
genesis and dark matter candidates [35]. If the baryon asymmetry is generated (or dark matter 
density is fixed) before the curvaton dominates the Universe, too large isocurvature perturbation 
will be produced. Thus, both baryon asymmetry and dark matter must be generated after the 
curvaton domination. In the axionic curvaton, the Hubble parameter at the curvaton domination, 
-ffdom) increases as rrio- and /. For rn-o- = ITeV - 10^ TeV and / ~ 10^''GeV, it ranges from 1 GeV 
to 100 TeV. There are several baryogenesis mechanisms which work at a Hubble parameter below 
Hdom- For instance, in the Affleck-Dine mechanism [36, 37], the baryon number is generated and 
fixed when the Hubble parameter is comparable to the soft mass of the flat direction in the MSSM. 
For the sfermion masses of order 10 - 100 TeV, it is possible that the AD field starts to oscillate 
after the curvaton domination. Since the mass of the AD field at large field value has rather large 
uncertainty, i^dom below TeV may be also allowed; for instance, this is the case if the potential of 
the AD field becomes flatter at large fields values. There are many dark matter candidates. Since 
the curvaton decays just before BBN for the case of our interest, there is an entropy dilution. One 
of the plausible dark matter candidates is the QCD axion, which starts to oscillate when the plasma 
temperature drops down to the QCD scale. There may be other ultralight axions which contribute 
to the dark matter density. Also, the thermal relic abundance of the WIMPs as well as WIMPs 
non-thermally produced by the curvaton decays are candidates for the dark matter. 

We have assumed that the curvaton is responsible for the observed density perturbation. From 
the minimalistic point of view, of course, the quantum fluctuation of the inflaton is the leading 
candidate. However, requiring both an extremely flat potential for sufficiently long inflation and the 
normalization of density perturbation may be too strong constraint on the inflation sector. If there 
are many other light scalars in nature, it might be more probable that there are two scalars, namely, 
the inflaton and the curvaton, responsible for the inflationary expansion and the origin of density 
perturbations, respectively. 

5 Conclusions 

In this paper we have studied non-Gaussianity of the density perturbation generated by the axionic 
curvaton, focusing on the case that the curvaton initially sits near the hilltop of the potential during 
inflation, and dominates the Universe before it decays. Interestingly, we have found that the non- 
Gaussianity parameter /nl is positive and gets enhanced up to 30 (or 40 for early reheating) in 
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the hilltop limit, even when the curvaton dominates the Universe. We have confirmed that this 
conclusion holds for = 0.94-0.99. It was also shown that in extreme cases where the axionic 
curvaton drives a secondary inflation, then the produced non-Gaussianity is typically /nl ^ 10 and 
is smaller than non-inflating cases. Note that, as long as the curvaton dominates the Universe, 
/nl cannot be larger than 30-40; this should be contrasted to other scenarios which can generate 
arbitrarily large non-Gaussianity, and some parameters must be tuned to realize /nl = 0(10). We 
have also pointed out that one of the plausible candidates for the axionic curvaton is the string axion 
with mass of order 10- 100 TeV and decay constant of 10^^"^^ GeV. If there are many axions in the 
Universe, one of them may be indeed responsible for the origin of the density perturbation. 
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A Inflating Curvatons 

In this appendix we consider the possibility that the curvaton drives a second inflationary stage 
before it starts to oscillate. After giving general discussions on density perturbations sourced by 
such inflating curvatons, we study the case for axionic curvatons. 

A.l Density Perturbations from Inflating Curvatons 

The case studied in this appendix is illustrated in Figure 2: The curvaton initially (i.e. during 
the first inflationary era) has negligibly tiny energy density compared to the total energy, however 
dominates the universe before it starts its oscillation. We suppose that this second inflationary 
period is not so long, and the CMB scale exits the horizon during the first inflation. The curvaton's 
field fluctuations obtained from the first inflation lead to slight difference in the lengths of the second 
inflationary periods among different patches of the universe, thus generate the density perturbations. 
We also note that the second inflation is not necessarily a slow-roll one, but may be a rapid- roll 
inflation [38, 39, 40, 41, 42], depending on the curvaton potential. 

Upon calculating the density perturbations using the (5AA-formalism, we assume that the second 
inflation lasts long enough (say, more than one e-fold) such that during this period the inflaton 
energy density becomes negligibly tiny and then the universe is well described as composed only of 
the curvaton. This assumption allows us to choose the final uniform energy hypersurface to possess 
energy density pf equal to or larger than that at the end of the second inflationary period, cf. 
Figure 2. In other words, we set the final hypersurface to be before the end of the curvaton-driven 
inflation, but late enough so that after which the inflaton can be ignored. 

^^Note that the curvaton field fluctuations at the CMB scale is generated during the fast inflation, thus the resulting 
density perturbations can be (nearly) scale-invariant even if the second inflation is a rapid-roll one. 
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Figure 2: Scliematic of tlie time variation of energy densities in an inflating curvaton scenario. 

The energy density of the inflaton (/> is considered to redshift as oc a^^ after the first infla- 
tion, and we first study the case where the inflaton decay happens after the curvaton domination. 
Moreover, the curvaton a is assumed to drive slow/rapid-roll inflation after dominating the universe. 
Hereafter we use the subscripts * to denote values when the CMB scale exits the horizon (during 
the first inflation), "end" for values at the end of inflation, "dom" for when the curvaton dominates 
the universe (i.e. /O^dom = Po-dom), and "/" at the final constant energy density hypersurface. Then 
in order to compute the resulting density perturbations, we would like to obtain the o"^,-dependence 
of the e-folding number from the end of inflation until the final surface {p^ is negligibly tiny during 
the first inflation, thus the curvaton has little effect on the expansion history before tend)- 

M = Na+Nb, (A.l) 

where 

Na= I Hdt, Mb= Hdt. (A.2) 

Here H = d/a, with an overdot denoting a time-derivative. 

We take V{a) to be the energy potential of the curvaton, which we assume to have no explicit 
time dependence. Then using p^dom — ^(o'dom)) one finds 

Ma = \ln^^\ln-^^ (A.3) 

o Pcfidom o V {adorn) 

thus 

dMg 1 ^'(o-dom) dadom /^^x 

da* 3 V{adom) da* 
where a prime denotes a derivative in terms of a. However, we will soon see that 6Ma only gives a 
minor contribution to the density perturbations. 

After the curvaton domination, for simplification, we ignore /O^ and describe the second inflation 
as a single-component slow/rapid-roll inflation. Then the inflationary dynamics is approximated 



^''This approximation is valid as long as the main contribution to SAf comes from the difference in the duration of 
the second inflation. 
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by (cf. appendix of [42] ) 

^M^H^c^V, cHa^-V', where c= ^ + VQ- ^g!?^ ^ = ^2^, (A.5) 

which are stable attractors under the condition 

Ml fV'V 

and for a nearly constant rj satisfying r] < 3/4. Here, note that ry is not bounded from below, and 
that c > 3/2. The familiar slow-roll approximations are recovered when |r/| <^ 1. To be precise, 
inflation can happen even for e > 1 given a large c (i.e. largely negative r/).^'' However in this 
appendix we limit our studies to curvaton potentials satisfying (A. 6) at cr = (Tdomj as the hilltop 
potentials which are discussed in the next section satisfy this condition. Then, since the curvaton 
field value is monotonically increasing or decreasing in terms of time, we can use a as a clock, 

rCTf TT 

Nb= —da. (A.7) 

Here, note that H/& is a function of cr, and since pf is a constant among different patches of the 
universe, so is Cf}^ Hence by partially differentiating both sides in terms of cr*, one obtains 

dMb cV 



— ""dom 



In order to compute daAoia/ da^,, we make use of the slow-roll approximation 2>Ha ~ — while 
t < tend- Then during tend < i < idom, for simplification we treat the universe as a matter dominated 
one and adopt \H(t ~ —V' (cf. (2.1), see also Footnote 16). Also using H / H"^ = —3/2 for 
tend < t < tdom, then one can check that 

/ 777^^ — / -ppr + (terms independent of cr*). (A. 9) 

Partially differentiating both sides by cr*, and using 3MpH^^^ ~ 2V{adom),^^ one obtains 

da^ ^ V {adorn) ( ) 

5(7* y'((j*) ' ^ ^ 

where we have dropped the contribution from the right hand side of (A. 9) from the condition (A. 6). 



Combining the above results, we can calculate the density perturbation spectrum: 

/ r/zv \ I n ... \^ 



^^The study on the stability of the rapid-roll attractor given in the appendix of [42] mainly considers c — 0{1), but 
one can easily extend their discussions to cases with c 3> 1. 

^^Considering pf ~ V{af), then for a potential that monotonically increases or decreases in terms of a during 
inflation, one sees that a/ is a constant. 

^^This may seem contradicting with the slow/rapid-roll approximation (A.5) 3MpH^ ~ V, but the numerical 
coefficient of V only affects a Mp{V' /V)'^ term which is dropped in the final expression (A. 10), thus we will not worry 
about it. 
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where (again using (A. 6)) 



The spectral index follows as 



taking the same form as for non-inflating curvatons (2.9). The non-Gaussianity parameter can also 
be calculated: 

6dal\daJ "2c>dom)\V ^(^dom) J n'^dom) J ' ^ ' ' 

where • • • denotes terms proportional to dc{adom) / dadom- One immediately sees that the first term 
in the { } parentheses is much smaller than unity from the condition (A. 6), while the second term 
can be larger than unity for rapid-roll inflation. 



_ c(adom)^(c^dom) -^2) 



H* 2V"iaA , 



In the above discussion, we have considered the inflaton to decay after the curvaton domination. 
Similar computations can be carried out also for the case where the inflaton decays between the 
first and second inflationary periods, by approximating the universe as matter dominated while 
tend < i < ireh (here the subscript "reh" denotes values H = T^, when the inflaton is assumed 
to suddenly decay), and then radiation dominated while treh < t < tdom- Further assuming the 
condition (A. 6), and also that the tilt of the curvaton potential at cjreh to be not much greater than 

f-t Cdomj i-6., 

\V'{<yr,Y.)\<\V'{adom)\, (A.15) 

then one obtains the same results (A. 12), (A. 13), and (A. 14). Whether the inflaton decays before 
or after the curvaton domination has little effect since the density perturbations are sourced mainly 
through different patches of the universe experiencing slightly longer /shorter periods of the second 
inflation. 



In summary, independently of whether the inflaton decays before/after the curvaton domination, 
under the condition (A. 6) (and also (A.15) for treh < ^dom); the linear perturbation sourced by 
inflating curvatons is of the form (A. 11) with (A. 12), the spectral index is (A. 13), and the non- 
linearity parameter is given by (A. 14). 

A. 2 Inflating at the Hifltop 

As an example, let us consider inflating curvatons with a hilltop potential 

V{o) = Vo-\m\a-a^)\ (A.16) 



where Vb, m, and ctq are constants. Given that the curvaton is located sufficiently close to the 
hilltop (To such that Vq ^ m?'{a — uo)^ and Vq ^> Mpm'^{cj — o-q)^ until the curvaton dominates the 
universe, then one finds 

§^^-^±4^— ^. (A.17) 
acr* brj ((To — a*) 
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and 

" J = 12 ("^ + , (A.18) 

where 

(A.19) 

Note especially that rj is (almost) a constant which is negative in this example. We also remark that 
we have dropped the contribution on /nl from the first term in the parentheses of (A. 14) which is 
clearly smaller than unity. 



A. 3 Axionic Curvatons 

In this final subsection, we look into axionic curvatons infiating at the hilltop of the potential (3.1). 
The discussions in Section A. 2 can be applied to this case by simply substituting 

Vb = 2A^ "i' = ^, cTo = /vr. (A.20) 

Then one can see that the negative ry parameter is determined merely by the symmetry breaking 
scale / as 

^=^-2^- (A.21) 

Non-Gaussianities from inflating axionic curvatons are in general smaller compared to non- 
inflating cases, which can be seen from (A.18): Large /nl requires a largely negative r], however a 
large |??| substantially accelerates the curvaton thus makes it challenging for the curvaton to drive 
inflation in the first place. We show this explicitly in Figure 3, which investigates the parameter 
space of axionic curvatons beyond the allowed window for non-inflating axionic curvatons discussed 
in Section 3. Here we fix the curvaton mass to mg- = A^// = 10^ GeV, and investigate large / values 
beyond the right blue edge (corresponding to the requirement that the curvaton be subdominant 
until it starts oscillating) in Figure 1. Using (A. 13), we fix the (first) inflation scale Hi^f from the 
spectral index Ug ~ 0.96 assuming constant H during inflation, and also fix the curvaton position at 
CMB scale horizon exit cr*// from ~ 2.4 x 10^^ using (A. 17). Furthermore, we set the number of 
e-foldings between the CMB scale horizon exit and the end of the (first) infiation as A/"* = 50, the cur- 
vaton decay rate F^ by (3.3) with /? = 1, and the infiaton decay rate F,^ small enough such that the 
infiaton decays after the curvaton domination. (The explicit value of F^ is irrelevant for the density 
perturbations, however whether the infiaton decays before/after the curvaton domination slightly 
affects the number of e- folds obtained in the second inflation.) The resulting non-Gaussianity /nl 
is plotted as a function of / in Figure 3, where the blue solid line denotes the analytic calculation 
(A.18) with (A.21). We have also numerically computed /nL) whose results are shown as blue dots 
in the flgures. One sees that the analytic and numerical results match well. In the right flgure, we 
also show the number of e-folds A4ec obtained in the second inflationary period driven by the axionic 
curvaton. Msec here is deflned as the e-folding number from the curvaton domination until when the 
curvaton starts oscillating, i.e. (2.2). When the second inflationary period is very short, the analytic 
estimations derived in this appendix are invalid, which sources the slight difference between the 
analytic and numerical computations of /nl at / « 10^^'^ GeV. When further increasing / beyond 
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the plotted regime, /nl becomes further suppressed while Msec rapidly increases, soon making the 
axionic curvaton responsible for driving most of the inflationary e-folds after the CMB scale horizon 
exit. In summary, non-Gaussianity in the region beyond the right blue edge in Figure 3 decreases for 
larger /, taking values smaller than ~ 10 in most of the region. Larger non-Gaussianity is generated 
when closer to the edge, i.e. when the second inflationary period is very short and the situation is 
close to the familiar non-inflating curvatons. 

Let us also note that the power spectrum (A. 17) is now written as 



r. 



1/2 

c 



1 



where 



3 + V9 - Uf] 



(A.22) 



For a sub-Planckian / (i.e. / < Mp), the prefactor k can only take values within 0.04 < k ^ 0.12. 
Therefore, once the initial position of the curvaton o"^,/7r/ is given, the inflationary scale -ff^, needs to 
be tuned to a rather narrow scale range in order for an inflating axionic curvaton to source the linear 
perturbation with an appropriate amplitude. This is in contrast to non-inflating axionic curvatons, 
which can work with a wide range of inflationary scales for each value of o'*/7r/ [10]. 
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Figure 3: Left: Non-Gaussianity from an inflating axionic curvaton with a fixed mass mo- = 10^ GeV, 
when increasing the decay constant / beyond the right blue edge of the allowed parameter window in 
Figure 1. Parameters other than m„ and / are fixed from requiring Vt; ^ 2.4 X 10-9 and ^ 0.96. 
Right: Number of e-folds in the second inflationary period driven by the axionic curvaton. 

We remark that the inflating curvaton was studied in Ref. [43] , however their results differ from 
ours. In particular, they claimed that the non-Gaussianity parameter /nl is negative and of order 
unity while we have shown that /nl from a hilltop curvaton is positive and can take values larger 
as well as smaller than order unity. Moreover we have confirmed our results for the case of axionic 
curvatons by numerical calculations as shown in Fig. 3, where one sees that /nl varies from 9 to 1 
as the e-folding number in the second inflation increases from 1 to 10. 



They assumed the curvaton mass to be much lighter than the Hubble parameter during inflation, but even in this 
case the analysis shown in Figure 3 does not change much since /nl from an axionic curvaton is set merely by the 
decay constant /. 
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